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A  conference  on  "Applied  Probability  Models  for  Complex  Stochastic 
Systems"  was  held  at  Jane  S.  McKimmon  Center,  North  Carolina  State  University, 
on  January  8-9,  1981.  The  conference  was  organized  by  Dr.  Shaler  Stidham, 

Jr.,  Professor  in  the  Department  of  Industrial  Engineering  and  Graduate 
Program  in  Operations  Research,  North  Carolina  State  University. 

The  purpose  of  the  conference  was  to  bring  together  several  prominent 
mathematicians  and  engineers  from  Europe  and  the  U.S.  with  common  interests 
in  applied  probability  models  for  complex  stochastic  systems,  to  interchange 
ideas  with  one  another  and  with  faculty,  students,  and  researchers  at  the 
universities  and  research  organizations  in  the  Research  Triangle  area. 

Primary  emphasis  was  placed  on  two  topics:  (i)  distribution-free 
analysis  of  complex  stochastic  systems;  and  (ii)  optimal  control  of  complex 
stochastic  systems.  This  emphasis  was  motivated  by  the  need  for  mathematical 
models  to  help  describe  and  control  such  complex  stochastic  systems  as  networks 
of  queues,  such  as  arise  in  computer  and  communication  systems  and  industrial 
job  shops,  and  operating  systems  composed  of  parts  subject  to  maintenance, 
failure,  and  replacement,  such  as  arise  in  aerospace  and  military  applications. 
(More  detailed  discussion  of  the  topics  covered  and  the  rationale  for  choosing 
them  can  be  found  in  the  proposal  for  the  conference,  which  is  Appendix  A 
to  this  report.) 

The  speakers  at  the  conference  were: 

Dr.  David  Burman,  Bel  1  Laboratories 

Prof.  Awi  Federgruen,  Columbia  University 

Prof.  Robert  Foley,  Virginia  Polytechnic  Institute 

Prof.  Peter  Franken,  Humboldt  University,  Berlin  (visiting  Virginia 
Polytechnic  Institute) 

Dr.  Oaniel  P.  Heyman,  Bell  Laboratories 

Prof.  Frank  Kelly,  Cambridge  University 

Prof.  S.  S.  Nair,  St.  Augustine's  College,  Raleigh 
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Prof.  N.  U.  Prabhu,  Cornell  University  /  Ar 

i  A/ 

/  .  '  >’ 

Prof.  Tomasz  Rolski,  Wroc7aw  University,  Poland  !  !'  ■  . , 

/ 

Prof.  Rolf  Schassberger ,  Technical  University,  Berlin  -  X) 

Dr.  Richard  F.  Serfozo,  Bell  Laboratories  >  Kl?.  ^ 

Dr.  Burton  Simon,  Bell  Laboratories  , .  ‘ 

/  *  7 
; 

Dr.  D.  R.  Smith,  Bell  Laboratories  - 

Prof.  Henk  Tijms,  Free  University,  Amsterdam  \  fit 

Prof.  Jean  Walrand,  Cornell  University 
Prof.  Jaap  Wessels,  Eindhoven  University,  The  Netherlands 
A  copy  of  the  program  for  the  meeting,  including  titles  of  the  talks, 
is  included  in  Appendix  B  to  this  report.  Detailed  abstracts  and,  in  some 
cases,  complete  papers  for  the  talks  are  included  in  Appendix  C. 

We  wish  to  express  our  gratitude  to  each  of  the  speakers  for  the  time  and 
effort  that  went  into  their  talks,  which  contributed  to  the  success  of  the 
meeting.  Thanks  are  also  due  to  Prof.  Ralph  Disney,  who  attended  the 
meeting,  and  to  the  other  organizers  of  the  TIMS  meeting  on  "Applied  Probability- 
Computer  Science:  The  Interface",  held  January  5-7,  1981,  at  Florida  Atlantic 
University,  Boca  Raton,  Florida.  Most  of  the  speakers  at  our  conference  also 
attended  that  meeting,  without  which  it  would  not  have  been  possible  to  bring 
together  so  many  distinguished  researchers  in  Raleigh.  We  are  especially 
grateful  to  Prof.  Disney  for  his  encouragement,  helpful  suggestions,  and 
cooperation,  particularly  in  making  it  possible  for  Prof.  Franken  to  attend. 

In  addition  to  the  speakers,  there  were  about  30-35  participants  from  the 
Research  Triangle  area  and  elsewhere.  These  included  the  following  persons: 

Prof.  John  W.  Bishir,  North  Carolina  State  University 

Prof.  Halim  Dogrusoz,  North  Carolina  State  University 

Prof.  S.  E.  Elmaghraby,  North  Carolina  State  University 

Prof.  George  S.  Fishman,  University  of  North  Carolina,  Chapel  Hill 

Prof.  Arnold  Greenland,  George  Mason  University 
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Or.  Robert  L.  Launer,  Army  Research  Office 

Prof.  Ross  Leadbetter,  University  of  North  Carolina,  Chapel  Hill 
Prof.  Louis  Moore,  University  of  North  Carolina,  Chapel  Hill 
Prof.  Arne  Nilsson,  North  Carolina  State  University 
Prof.  Henry  L.  W.  Nuttle,  North  Carolina  State  University 
Prof.  Walter  L.  Smith,  University  of  North  Carolina,  Chapel  Hill 
Prof.  M.  Venkatesen,  Case  Western  Reserve  University 
Dr.  Gerald  Andersen,  Army  Research  and  Development  Command 
We  would  like  to  thank  the  Army  Research  Office  for  providing  the  support 
which  made  this  conference  possible.  Special  thanks  should  go  to  Dr.  Robert 
L.  Launer  of  ARO,  for  his  personal  support  and  hi s efforts  to  expedite  the 
proposal  on  unusually  short  notice. 

Finally,  our  sincere  thanks  to  Prof.  Salah  Elmaghraby ,  Director  of  the 
Graduate  Program  in  Operations  Research,  for  his  encouragement  which  led 
to  the  idea  of  holding  the  conference,  for  his  unselfish  comnitment  of  time 
and  advice  to  help  make  it  a  success,  and  the  financial  support  provided  by 
the  Operations  Research  Program. 


APPENDIX  A 


PROPOSAL  FOR  CONFERENCE 


A  Proposal  for  a  Conference  on 


APPLIED  PROBABILITY  MODELS  FOR 
COMPLEX  STOCHASTIC  SYSTEMS 


Subject  of  Proposal 

A  conference  on  "Applied  Probability  Models  for  Complex  Stochastic 
Systems,"  to  be  held  at  North  Carolina  State  University,  January  8-10,  1981, 
with  5-10  invited  speakers  (prominent  researchers  from  Europe  and  the 
United  States)  and  10-20  additional  participants  (faculty  and  students  from 
the  Research  Triangle  universities  and  applied  researchers  from  Research 
Triangle  Park  and  local  industry). 

Conference  Theme  and  Topics 

The  subject  of  the  conference  will  be  applied  probability  models  for 
complex  stochastic  systems,  particularly  multi-dimensional  systems  of 
interrelated  components.  Examples  are  networks  of  queues,  such  as  arise 
in  computer  and  communication  systems  and  industrial  jobs  shops,  and 
operating  systems  composed  of  parts  subject  to  maintenance,  failure,  and 
repair  or  replacement,  such  as  arise  in  aerospace  and  military  applications. 
These  topics  are  cu  rrently  of  great  interest  both  to  applied  probabilists ,  who 
are  developing  descriptive  and  normative  models  for  such  systems,  and  to  users 
and  potential  users  of  these  models  it  industry  and  the  military. 

Rather  than  consider  a  broad  range  of  possible  models,  we  intend  to 
place  primary  emphasis  on  two  topics: 

(i)  dist r ibut ion-f ree  analysis  of  complex  stochastic  systems;  and 

(ii)  optimal  control  of  complex  stochastic  systems. 
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Under  the  first  topic,  we  shall  focus  attention  on  such  subjects  as 
conservation  relations,  the  insensitivity  phenomenon,  bounds,  approximations, 
and  monotonicity  of  systems.  Investigation  of  these  subjects  will  be  based 
primarily  on  the  theory  of  regenerative  processes  and  the  more  general 
random  marked  point  processes  (RMPP) ,  sample-path  (operational)  analysis, 
and  stochastic  ordering.  Under  the  second  topic,  we  shall  consider  the 
structure  of  optimal  control  policies  and  efficient  computational  algorithms, 
based  primarily  on  the  theory  of  Markov  decision  processes  (MDP)  and  the 
techniques  of  dynamic  programming. 

Rationale  for  Conference 

There  are  several  reasons  for  proposing  a  conference  on  these  topics  at 
this  time. 

First,  except  for  simple  exponential  models  (such  as  the  Jackson 
network),  multi-dimensional  models  in  queueing  and  reliability  theory 
are  notoriously  difficult  to  solve  by  conventional  techniques,  such  as 
Markov-process  theory.  The  theory  of  RMPP  has  provided  a  tool  for  studying 
such  systems,  without  the  requirement  for  exponential,  or  even  independent, 
interarrival,  service,  failure,  or  repair  time  distributions.  Together 
with  the  concept  of  Palm  probabilities,  this  theory  has  facilitated  a 
unified  approach  in  a  distribution-free  setting  to  such  problems  as  the 
relation  between  customer  and  time  averages  in  queues,  the  relation  between 
a  continuous-parameter  stochastic  process  and  an  imbedded  discrete-parameter 
process,  and  the  insensitivity  of  equilibrium  distributions  in  queueing 
networks.  Parallel, but  so  far  more  limited  results,  have  been  achieved  using 
sample-path  analysis.  Using  some  of  these  results,  together  with  the  concept 
of  stochastic  ordering,  researchers  in  applied  probability  have  been  able  to 
approximate  complex  stochastic  systems  by  simple  systems,  with  tight  upper  and 
lower  bounds  for  such  measures  as  expected  queue  lengths  and  waiting  times. 
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The  theory  of  MDP,  while  not  extending  the  range  of  models  that  can  be 
solved  analytically,  has  focussed  attention  on  optimization  of  stochastic 
systems,  often  revealing  the  structure  of  optimal  control  policies  and  hence 
indicating  the  direction  in  which  descriptive  modelling  should  proceed.  In 
addition,  recent  computational  developments  in  dynamic  programming  theory 
have  increased  the  potential  scope  for  efficient  numerical  computation  of 
optimal  control  policies,  even  for  multi-dimensional  systems. 

A  second  reason  for  choosing  these  topics  is  that  research  in  each  has 
been  developing  in  recent  years  at  a  rapid  pace  in  Europe  (primarily  West  and 

East  Germany  and  the  Metherlands) .  As  a  consequence,  applied  probabi lists 
in  the  U.S.  are  perhaps  not  as  aware  of  the  latest  developments  as  they 
might  otherwise  be. 

Finally,  at  the  time  of  the  proposed  conference,  several  researchers 
with  interests  in  these  topics  from  Europe  and  the  U.S.  will  be  returning 
from  the  conference  on  "Applied  Probability-Computer  Science,  The  Interface," 
sponsored  by  the  ORSA/TIMS  Special  Interest  Group/College  on  Applied 
Probability,  to  be  held  at  Florida  Atlantic  University,  Boca  Raton,  Florida, 
January  5-7,  1981.  It  will  be  convenient  for  some  of  them  to  stop  in  Raleigh 
after  that  conference.  (It  was  the  fact  that  several  of  the  proposed 
participants  contacted  us  and  indicated  a  desire  to  visit  North  Carolina 
State  University  on  the  way  home  from  the  Florida  conference  that  originally 
provided  the  impetus  for  organizing  a  conference  in  Raleigh.) 

While  the  topics  proposed  for  this  conference  will  also  form  part  of 
the  agenda  at  the  Florida  meeting,  the  two  conferences  should  be  viewed  is 
complementary  rather  than  competitive  in  scope  and  purpose.  The  theme  of  the 
Florida  meeting  -  the  application  of  appLied  probability  models  to  computer 
systems  -  will  be  one  of  several  potential  applications  of  the  theoretical 
topics  in  the  proposed  conference.  On  the  other  hand,  with  its  focus  on 
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distribution-tree  analysis  and  control  models  for  complex  systems,  the  proposed 
conference  promises  to  be  more  narrow  in  scope.  Moreover,  we  expect  to  have 
a  much  smaller  number  of  speakers  (5-10  vs.  40-50)  and  a  more  homogeneous 
group  of  participants,  with  more  opportunities  to  pursue  the  selected 
theoretical  topics  in  depth. 

Benefits  of  Conference 

The  conference  will  take  place  over  a  period  of  2  1/2  -  3  days  and  will 
be  organized  around  5-10  invited  papers,  with  ample  time  for  discussion  of  the 

papers  and  informal  communication  among  the  participants.  We  hope  for 
two  main  benefits  to  come  from  the  conference.  The  first  will  be  the 
mutual  stimulation  from  bringing  together  a  small  group  of  applied 
probabilists  of  similar  interests,  each  of  whom  has  contributed  significantly 
to  the  development  of  models  for  complex  stochastic  systems.  The  second 
will  be  the  dissemination  of  knowledge  from  the  invited  speakers  to  the  applied 
probability  community  in  the  Research  Triangle  area,  including  faculty  and 
students  at  the  area  universities  as  well  as  applied  researchers  from 
Research  Triangle  Park  and  local  industry.  There  is  an  active  group  of 
faculty  and  students  at  North  Carolina  State  University  who  have  been 
contributing  substantially  in  recent  years  to  research  on  applied  probability, 
particularly  optimal  control  of  queueing  systems  and  modelling  and  numerical 
analysis  of  computer  and  communication  systems.  it  is  hoped  that  the  conference 
will  help  bring  recognition  and  added  stimulus  to  the  efforts  of  this  group. 

Possible  Participants 

A  tentative  list  of  invited  speakers  follows.  (Those  marked  bv  a  * 
have  been  contacted  and  have  indicated  a  desire  to  come.) 


*Ralph  Disney,  Virginia  Polytechnic  Institute  and  State  University 
*Peter  Franken,  Humboldt  University,  Berlin 


Daniel  P.  Hevman,  Bell  Telephone  Laboratories 
Frank  Kelly,  Cambridge  University 
*Dieter  Konig,  Bergakademie  Freiberg 

*Austin  Lemoine,  Systems  Control,  Inc.,  and  Stanford  University 
*Teunis  Ott,  Bell  Telephone  Laboratories 
*N.  U.  Prabhu,  Cornell  University 
*Tomasz  Rolski,  Wroclaw  University 
Rolf  Schassberger ,  Technical  University,  Berlin 
*Richard  Serfozo,  Bell  Telephone  Laboratories 
*Henk  Tijms,  Free  University,  Amsterdam 
*Jaap  Wessels,  Eindhoven  University,  The  Netherlands 
Ward  Whitt,  Bell  Telephone  Laboratories 

The  following  is  a  partial  list  of  researchers  from  the  Research 
Triangle  area  who  might  also  participate. 

B.  B.  Bhattacharyya,  North  Carolina  State  University 

John  W.  Bisher,  North  Carolina  State  University 

W.  W.  Chou,  North  Carolina  State  University 

Halim  Dogrusoz,  North  Carolina  State  University 

George  Fishman,  University  of  North  Carolina,  Chapel  Hill 

R.  Leadbetter,  University  of  North  Carolina,  Chapel  Hill 

Arne  Nilsson,  North  Carolina  State  University 

Walter  Smith,  University  of  Nortii  Carolina,  Chapel  Hill 

William  R.  Stewart,  North  Carolina  State  University 


Other  interested  participants  will  also  he  welcomed. 
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BUDGET 


January  8-10,  L  98  L 


Requested  From  NCSU  and 

from  ARO  Other  Sources 


*1) 

Invited  Speakers:  travel  and  living 
expenses  @  $400/person 

$  2400 

$ 

2800 

2) 

Preparation  of  Papers 

1000 

— 

3) 

Local  Participants:  travel  and  meal 
expenses  for  10  persons 

300 

-- 

4) 

Facilities:  room  charges,  etc., 

at  McKimmon  Center 

2 1)0 

-- 

3) 

Miscellaneous:  postage,  telephone,  etc. 

2  30 

-- 

*6 ) 

Travel  Expenses:  for  Dr.  F.  P.  Kelly 
Round  trip  airfare  England-U.S. 

(300 

— 

TOTAL 

S  4930 

S 

28on 

'■'AIL  the  invited  participants,  with  the  exception  of  Dr.  Kelly,  expect  to 
receive  support  from  other  sources  to  attend  the  Florida  TIMS/ORSA  meeting 
("Applied  Probability  -  Computer  Science,  The  Interface,"  January  3-7,  1931)  . 
Consequently  we  are  requesting  only  incremental  travel  expenses  to  and 
from  Boca  Raton,  Florida,  for  those  invited  speakers  who  need  such  support. 

In  the  case  of  Dr.  Kelly,  however,  we  are  requesting  support  for  round-trip 
travel  from  England  to  the  U.S.  as  well,  because  we  have  been  informed  by 
Dr.  Teunis  Ott,  Organizing  Chairman  of  the  FLorida  meeting,  that  without  such 
support  he  will  not  be  able  to  attend  either  meeting. 

NOTE:  No  money  from  ARO  will  be  used  to  invite  speakers  from  Eastern  Block 

countr  ies . 


Civil  Rights  Compliance  Statement 

"No  person  will,  on  the  grounds  of  race,  color  or  national  origin  be  excluded 
from  participation,  be  denied  the  benefits  of,  or  be  subjected  t  o  discrimin.it  ion 
under  this  program." 

A  f  f  i  rma t ive  Action  Statemen t 

"NCSU  does  not  discriminate  on  the  basis  of  sex,  race  or  handicap  and  is 
an  affirmative  act ion/eoua  1  employment  institution." 
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REVISED  PROGRAM 


Conference  on  Applied  Probability  Models 
for  Complex  Stochastic  Systems 

January  8-9,  1981 


McKimmon  Conference  Center 
North  Carolina  State  University 
Raleigh,  NC 


Schedule  (subject  to  change) 

(Shuttle  bus  will  leave  Mission  Valley  Inn  at  8:30  each  morning,  return 
from  McKimmon  Center  at  5:30  each  evening.) 


Thursday,  January  8 

8:45  -  9:00  Registration 


9:00  -  9:15 
9:15  -  10:00 
10:00  -  10:45 
10:45  -  11:00 
11:00  -  11:45 
11:45  -  12:30 
12:30  -  2:00 
2:00  -  2:45 

2:45  -  3:30 

3:30  -  3:45 


Opening  remarks 

R.F.  Serfozo,  "Compound  Point  Processes  and  Rare  Events" 
Rolf  Schassberger,  "Networks  of  Queues  in  Discrete  Time" 
Coffee  break 

Frank  Kelly,  "Networks  with  Partial  Balance" 

Jean  Walrand,  "Filtering  Theory  and  Networks  of  Queues" 
Lunch* 

David  Burman,  and  D.R.  Smith,  "A  Light-Traffic  Theorem 
for  Multi -Server  Queues" 

N.U.  Prabhu,  "Comparison  of  Random  Walks,  with  Applications 
to  Queues" 

Coffee  break 


3:45  -  4:30  Robert  Foley,  "Compartmental  Models  and  Marked  Point 
Processes" 
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4:30  -  5:15  S.  S.  Nair,  "A  Carrier  Shuttling  between  Two  Terminals" 
6:30  Cocktails  and  dinner  (Balentine's  Restaurant,  Cameron  Village)** 


Friday,  January  9 
9:00  -  9:15 
9:15  -  10:00 

10:00  -  10:45 

10:45  -  11:00 
11:00  -  11:45 

11:45  -  12:30 
12:30  -  2:00 
2:00  -  2:45 


Announcements 

Henk  Tijms,  "Stochastic  Control  for  Specially  Structured 
Optimization  Problems" 

Awi  Federgruen,  "Approximations  of  Dynamic,  Multi -Location 
Production  and  Inventory  Problems" 

Coffee  break 

Jaap  Wessels,  "Using  the  Problem  Structure  in  Large 
Markov  Decision  Processes" 

Peter  Franken,  Topic  to  be  announced. 

Lunch* 

Daniel  Heyman,  "Conservation  Laws  for  Queues" 


2:45  -  3:30 


3:30  -  3:45 


Tomasz  Rolski,  "Queues  with  Non-Stationary  Input  Stream: 
Ross's  Conjecture" 

Coffee  break 


3:45  -  4:30 


Burton  Simon,  "Priority  Queues  with  Feedback" 


*  Included  in  registration  fee  of  $20.00  (or  $10.00  for  one  day).  Spaces 
are  limited.  Registration  fee  for  conference  only  (no  lunch):  $4.00 
per  day  (no  charge  for  NCSU  students). 

**  Dinner  reservations:  $10. 00/person.  Spaces  are  limited. 


For  further  information  contact: 

Professor  Shaler  Stidham 
North  Carolina  State  University 
Department  of  Industrial  Engineering 
Box  5511,  Raleigh,  NC  27650 
(919)733-2362 
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Abstract  of 

Compound  Point  Processes  and  Rare  Events 

by 

Richard  Serfozo 

Bell  Laboratories,  Holmdel  N.J. 

It  is  well  known  that  the  binomial  distribution  for 
rare  events  converges  to  (and  hence  can  be  approximated 
by)  a  Poisson  distribution.  I  will  show  how  this  conver¬ 
gence  result  for  distributions  of  rare  events  extends  to 
Doln.t  processes.  First,  a  "compound"  point  process  for  events 
with  attendant  random  variables  will  be  defined  usinsj 
Laplace  functionals.  Then  several  theorems  will  be  present¬ 
ed  that  characterize  the  convergence  of  these  compound 
processes  as  their  defining  parameters  tend  to  values 
indicative  of  rare  events.  Applications  to  partitioning 
and  merarlna;  of  streams  of  customers  in  networks  will  also 
be  discussed. 
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Networks  of  Queues  in  Discrete  Time 
by 

Rolf  Schassberger 
Technical  University 
Berlin 

ABSTRACT 

The  networks  of  queues  described  and  analyzed  in  this  talk  consist  of 

so-called  doubly  stochastic  servers.  When  in  isolation  such  a  server  operates 

in  discrete  time  as  follows:  the  time  is  sliced  into  quanta  of  equal  length; 

during  each  quantum  at  most  one  job  may  arrive  at  the  server,  the  probability 

for  an  arrival  being  > ,  A  e  (0,1);  an  arrival  is  of  type  r,  r  ,  R,  with 

probability  a  ,  where  R  is  a  countable  set  and  l  ct  =  1 ;  an  arrival  of  type 
r  r.:R  r 

r  carries  a  request  of  i  quanta  of  service  time,  ,.r  =  1,2,...;  the  server 
hands  out  one  quantum  of  service  to  each  job  immediately  upon  arrival  and 
thereafter  operates  under  a  so-called  doubly  stochastic  discipline;  a  complete 
description  of  these  disciplines  being  fairly  involved,  it  may  be  sufficient 
here  to  state  that  they  include  round-robin,  last-in-fi rst-served-preemptive- 
resume,  and  random-service  scheduling,  all  with  the  additional  feature  of 
queue-size  dependent  service  capacity;  upon  completion  of  its  request  a  job 
departs  from  the  server. 

It  is  clear  from  this  description  that  the  doubly  stochastic  server  is  a 
rather  important  model.  For  instance,  under  certain  specifications  of  the 
parameters,  it  becomes  essentially  a  discretized  M/G/l  queue  under  round-robin 
discipline.  The  state  of  the  doubly  stochastic  server,  suitably  defined  and 
observed  at  the  successive  beginnings  of  the  time  slices,  fluctuates  as  a 
discrete-state  Markov  chain.  The  steady-state  distribution  of  this  chain  is 


found  and  turns  out  to  be  of  a  simple  product-form.  A  similar  result  for  the 

M/G/l  round-robin  queue  is  not  available.  The  average  response  time  of  a  job 

of  type  r  is  found  as  well  and  turns  out  to  be  linear  in  . 

r 

Now  it  is  shown  how  doubly  stochastic  servers  can  be  linked  together  to 
form  networks  whose  modelling  power  is  similar  to  the  well-known  networks 
studied  by  Baskett,  Chandy,  Muntz,  and  Palacios,  and  by  Kelly.  Furthermore 
it  is  shown,  for  these  networks,  that  a  steady-state  law  of  product-form 
and  a  law  for  the  average  response  time  of  a  job  along  a  given  path  can  be 
easily  derived.  It  is  concluded  that  the  study  of  these  networks  should  be 
at  least  as  useful  as  that  of  their  well-known  counterparts . 
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Partial  balance 


by 

F.P.  Kelly 

Koenig,  Matthes  and  Nawrotski  (1967)  have  presented 
a  detailed  analysis  of  the  interrelationship  between 
insensitivity  and  partial  balance  within  a  framework  later 
christened  by  Schassberger  (1977)  a  Generalized  semi-Markov 
scheme.  Here  we  show  that  the  concept  of  partial  balance 
can  be  defined  for  an  unstructured  Markov  process  and  is 
then  equivalent  to  various  other  properties  of  the  process 
defined  in  terms  of  truncation,  time  reversal,  speed  changes 
and  observation  at  entrance  and  exit  times  (Kelly,  1979) . 
Finally,  we  show  how  the  relationship  partial  balance  implies 
between  the  invariant  measure  of  a  truncated  process  and  the 
invariant  measures  of  jump  chains  observed  at  entrance  and 
exit  times  can  be  used  as  the  basis  for  a  sample  path  con¬ 
struction  of  a  stochastic  process  exhibiting  insensitivity. 
This  approach  to  insensitivity  avoids  the  requirement  that 
the  successive  lifetimes  involved  in  the  definition  of  the 
process  be  generated  by  a  stationary  mechanism. 
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Filtering  Theory  and  Networks  of  Queues 
Jean  Walrand,  Cornell  University 

Abstract 

The  basic  idea  behind  the  filtering  approach  to  the  analysis  of  networks 
of  queues  is  to  "calculate"  the  conditional  distribution  ^  of  the  state  of 
the  network  given  the  past  (N$  ,0<^s£t)  of  an  observed  set  of  flows.  The 
filtering  theory  leads  to  the  system  of  (nonlinear)  differential  equations 
with  jumps  solved  by  A  related  system  of  linear  equations  solved  by  the 
un-normal ized  version  of  ^  can  also  be  derived. 

The  filtering  approach  provides  a  systematic  way  of  making  rigorous  some 
otherwise  delicate  limiting  arguments.  Some  extensions  of  the  arrival  theorem 
in  networks  of  quasi-reversible  node  constitute  an  example  of  those  applica¬ 
tions.  The  method  also  leads  to  useful  necessary  conditions  for  flows  in 
a  network  to  be  Poisson. 

Using  those  techniques  we  could  derive  the  characterization  of  Poisson 
flows  in  single  class  open  networks  of  quasi-reversible  nodes.  The  result 
is  the  intuitive  one:  A  flow  on  a  given  link  is  Poisson  if  and  only  if  the 
link  is  not  part  of  a  loop. 

Similar  ideas  can  also  be  used  to  obtain  the  conditions  under  which  two 
given  flows  in  a  network  have  the  same  law.  This  is  a  more  difficult  problem 
than  the  verification  of  the  Poisson  character  of  a  flow.  Only  a  few  non¬ 
trivial  examples  have  been  analyzed  so  far. 

Some  surDrising  equivalence  results  for  the  interarrival  time  distribu¬ 
tions  of  flows  in  networks  can  be  derived  by  combining  methods  from  filtering 
theory  and  first  step  analysis  in  Markov  chains. 

The  talk  concludes  with  a  remark  on  a  probabilistic  argument  for  the 
product  form. 


A  LIGHT-TRAFFIC  THEOREM 
FOE  MULTI-SERVER  QUEUES 


by 


David  Y.  Burman 
and 

Donald  R.  Smith 


Bell  Laboratories 
Holmdel,  New  Jersey  07733 


ABSTRACT 


Several  approximations  for  the  expected  delay  in  an  M/G/c 
queue  depend  on  both  its  light-  and  heavy-traffic  behavior. 
Although  the  required  heavy-traffic  result  has  been  proved, 
the  light-traffic  result  has  only  been  conjectured  by  Boxma 
Cohen,  and  Huffels.  We  prove  this  conjecture  when  the 
service  is  of  phase-type;  intuitively,  any  M/G/c  queue  car. 
be  arbitrarily  closely  approximated  by  such  a  system.  In 
particular,  as  the  traffic  goes  to  zero,  we  shew  that  the 
probability  of  delay  depends  only  or.  the  mean  of  the  servi- 
time  distributions  and  that  the  delay  (when  positive) 
converges  in  distribution  to  the  minimum  of  c  independent 
equilibrium-excess  service-times.  This  result  justifies  an 
efficient  computational  approach  to  obtain  numerical  result 
for  M/G/c  queues  and  provides  useful  methodology  for  the 
approximation  of  other  complicated  stochastic  systems. 
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Compartmental  Models  and  Marked  Point  Processes 

by 

R.  D.  Foley 

Industrial  Engineering  and  Operations  Research  Department 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  VA  24061 

A  compartmental  model  can  be  thought  of  as  a  system  of  m  locations  or 
compartments.  Particles  may  enter  the  system,  migrate  from  one  compartment 
to  another,  or  leave  the  system.  The  number  of  particles  in  each 
compartment  is  the  usual  quantity  of  interest.  For  example,  consider  the 
United  States  as  the  system,  each  state  a  compartment,  and  each  person  a 
particle.  As  people  are  bom,  die,  emigrate,  immigrate,  or  migrate  from 
one  state  to  another,  the  populations  of  the  states  fluctuate. 

Let  0  _<  _<  •  •  •  denote  the  sequence  of  arrival  times  of  particles 

to  the  system.  Associated  with  each  particle  is  a  random  element,  X^,  which 
describes  the  behavior  of  the  nth  particle  in  the  system.  X^  contains  the 
route  of  the  particle  and  the  length  of  time  the  particle  visits  each 
compartment  in  its  route.  Under  suitable  conditions  on  the  marked  point 
process  {(X^.T^):  n  >  0},  we  determine  the  joint  distribution  of  the  number 
of  particles  in  each  compartment  at  any  time  t. 

The  marked  point  process  {(X^.T^):  n  >  0}  is  required  to  satisfy  the 
following  two  assumptions: 

a.  The  number  of  arrivals  of  particles  in  the  interval 
[0,t]  forms  a  Poisson  process,  either  stationary  or 
non-s tat ionary; 

b.  The  behavior  of  any  particle  is  conditionally  independent 
of  the  behavior  of  the  other  particles  given  the  arrival 


time  of  the  particle.  That  is, 
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P{Xn  c  b|X1,T1,X2,T2, 


.X 


, , T  T  }  =  P{X 
-1  n-1  n  n 


e  B  T  }. 
'  n 


Define  Q^(t)  to  be  the  number  of  particles  in  compartment  i  at  time 
t,  D(t)  to  be  the  number  of  departures  from  the  system  in  [0,t],  and 
K_(t)  to  be  the  number  of  particles  that  migrate  from  compartment  i  to 
compartment  j  during  [0,t].  Under  the  assumptions  a  and  b  we  have  shown 
the  following  results: 

1.  Q, (t) ,Q~ (t) , • • • ,Q  (t)  are  independent  Poisson  random  variables; 

2.  {D ( t ) :  t  >  0}  is  a  Poisson  process  either  stationary  or  non¬ 
stationary; 


3.  Q^(t),Q2(t) ,  •  •  •  ,Q^(t)  is  independent  of  {D(u):  u  <_  t}; 

4.  (K„(t):  t  >  0}  is  a  Poisson  process  iff  almost  surely 
each  particle  migrates  from  i  to  j  at  most  once; 

5.  (K_(t):  t  >  0]  and  K^£(t):  t  .>  0}  are  independent 
Poisson  processes  iff  almost  surely  no  particle  migrates 
from  i  to  j  and  also  from  k  to  £,. 

Result  1  has  been  proven  for  tandem  compartmental  systems  in  Foley  (1981). 
Jackson  &  Aspden  (1980)  has  proven  result  1  for  tandem  compartmental  systems 
with  a  stationary  Poisson  arrival  process  and  exponential  compartmental  visit 
times.  Result  1  has  been  developed  for  a  particular  two  compartment  model  in 
Purdue  (1975).  Result  1  has  been  proven  for  a  single  compartment  model 
satisfying  assumptions  a  and  b  in  a  nice  paper  by  Brown  &  Ross  (1969).  Brown 
and  Ross  also  discuss  group  arrivals. 

Result  2  has  been  developed  for  a  single  compartment  system  with  constant 
visit  times  in  Ramakrishnan  (1980)  and  for  exponentially  distributed  visit 


times  in  Kambo  and  Bhalaik  (1979).  Result  2  is  related  to  Mirasol's  (1963) 


CIO 


result  that  the  departure  process  from  the  M/Gl/^  queue  is  a  stationary 
Poisson  departure  process.  In  Foley  (1980),  examples  are  given  of  queues 
with  non-stationary  Poisson  arrival  processes  with  stationary  Poisson 
departure  processes. 

Result  3  is  related  to  the  work  of  Kelly  (1979)  and  Melamed  (1979). 

The  property  in  result  3  appears  as  a  condition  in  both  of  these  papers 
in  order  to  obtain  the  analogue  of  result  2.  Result  3  is  quite  closely 
related  to  the  quasi-reversibility  condition  of  Kelly  (1979).  However 
it  is  not  exactly  the  same,  since  the  process  may  not  be  stationary.  The 
systems  studied  in  Foley  (1980)  possess  the  property  in  result  3  but  are 
non-stationary . 

Results  4  and  5  are  the  analogues  of  the  results  obtained  by  Melamed 
(1979)  and  Walrand  and  Varaiya  (1978)  for  Jackson  queueing  networks. 

The  results  obtained  for  the  compartmental  models  satisfying  conditions 
a  and  b  are  strikingly  similar  to  tne  recent  results  obtained  on  Jackson 
queueing  networks  and  variants  thereof  by  Kelly  (1979),  Melamed  (1979)  and 
Walrand  and  Varaiya  (1978). 
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i  tochasti^  coni:  1  fox  special  ly  stx  uctux  -  i  opt  j  mi  zatioii  proiiert.s 


'■■■■  .S'.  rSv .  In  the  lari  filter-:.  year:-  th«-i  «•  has  b<-on  a  cor.:.:  derable  interest  in 
stiiy  of  optimal  design  u:.<;  control  oi  production  at.  i  queu-. :  ng  system:  ,  cf. 
the  bibliography  by  Cra: ill  et  al .  (1977)  and  the  survey  papers  by  Prabhu 
-l.u  5 x - a;.e  Sooel  (ls7u).  On  the  one  hand  the  l.teiauii't  Ceam 
with  steady-state  analysis  of  an  intuitively  reasonable  control  rule  having 
a  simple  form  and  on  the  other  hand  a  large  number  cf  papers  are  cor.cerr.ec 
with  verifying  the-  optimality  of  such  a  simple  control  rule  among  a  la:  ger 
class  cf  control  rules.  Hovr-vcx,  so  far  little  attention  ha?  lee:,  paid  to 


Lem:  . 


mi -Mari-  o v  d c  i r.ion : 


'•  •-a-1  arV-  '  -1- 

Tii  .i  r. 


which  appear-o-I  t..  b<  quite  tesfull  for-  th'-  a;.;  1  i  cat  ions-  cons i -ier-t 
appro,  teb  cord  in--:;  policy-:  t  era  tier:  and  embedding  technique-?  in  such  a  way  that 
we  c-n.y  to  perform  c-ilcul a*  ions  fox-  «  ;V». r.um'.  er  cf  states  without 


states  in  any  iteration  and  gen 
the  desired  simple  form.  In  eu^. 


.  coals  only  with  a  finite  eel-- ice  :  set  oi 
rates  a  sequence-  of  improved  policies  having 
..  of  the  uppl icat ions  oonsiuexeu  ihe  uli.it- ns  ion 
of  the  embedded  sets  of  states  is  considerably  smaller  than  that  of  any 
appropriately  chosen  finite  state  space  approximation  and  this  dimensionality 
reduction  is  important  in  viev;  of  computations. 

After  having  discussed  th*-  embedding  approach  tc  compute  the  avei  ago  cost 
and  the  relative  cert  fui.ct we  pres-.-nt  a  numnex  oi  applic.it  ions .  The  first 
two  application:  do  a  i  wit:,  th*--  c-pt  inn  1  coi.no!  oi  queu.-i  a:..:  .  .  ncc*r 
queue  wit!,  variable  nu  vir-  rate  and.  an  f-l/M/n  queue  ir.  whl*  :.  tn*.  i.u 
operating  rex  vie-  channels  cm  !-•  vui  b-d.  !!<  :<i  w.  discu  s  r ,-m*  appl 
inventory  a:  :  pa  ..cue- :..  W*  :  ix  :  t  cons  id*  r  a  pi  iuc‘ 1c  ;  sy-  t  in 
production  c  curs  an:  : inv-ntoiv  c.<:  i  *.-  cm  trolled  iv  u;.; 
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production  facility.  An  efficient  algorithm  i  giv»-n  n  c  '-mp  u*e 
(2,s)  p.  l;.y.  V«>  finally  dic-our.  a  mult  -  i  tv.  •  y  '  -.-r.  in  which 
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■und  1’oi  r.sc'n  demands  an  ! 
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Approximations  of  dynamic,  multilocat i on 
production  and  inventory  problems 


by 

Awi  Federgruen 
Graduate  School  of  Business 
Columbia  University 
New  York,  NY  10027 

and 

Paul  Zipkin 

Graduate  School  of  Business 
Columbia  University 
New  York,  NY  10027 

ABSTRACT 

Consider  a  central  depot  (or  plant)  which  supplies  several  locations  exper¬ 
iencing  random  demands.  Orders  are  placed  (or  production  is  initiated)  periodic¬ 
ally  by  the  depot.  The  order  arrives  after  a  fixed  leadtime,  and  is  then  allocated 
among  the  several  locations.  The  allocations  are  finally  received  at  the  demand 
points  after  another  lag. 

This  system  gives  rise  to  a  dynamic  program  with  a  state  space  of  very  large 
dimension.  We  show  that  this  model  can  be  systematically  approximated  by  a 
single-location  inventory  problem,  employing  aggregation  techniques.  All  the 
qualitative  and  quantitative  results  for  such  problems  can  then  be  applied. 

KEYWORDS:  INVENTORY  6  PRODUCTION; 

APPROXIMATIONS, 


STOCHASTIC  MODELS 


Using  the  problem  structure  in  large  Markov  decision 


processes . 
by 

Jaap  Wessels 
Eindhoven 

This  talk  is  based  on  the  paper  "Some  notes  on  iterative  optimization  of 
structured  Markov  decision  processes  with  discounted  reward"  by  M.  Hendrikx, 

J.  van  Nunen  and  J.  We=';als  (memorandum  COSOR  80-20,  Eindhoven  University 
of  Technology,  Department  of  Mathematics,  November  1980) . 

The  paper  contains  a  comparison  of  solution  techniques  for  Markov  decision 
processes  with  respect  to  the  total  reward  criterion.  It  is  illustrated  by  4 
examples  (2  from  literature  and  2  from  our  own  experience) ,  that  the  effect  of 
a  number  of  improvements  of  iterative  methods,  which  are  advocated  in  the 
literature,  is  limited  in  some  realistic  situations. 

Numerical  evidence  is  provided  to  show  that  exploiting  the  structure  of  the 
problem  under  consideration  often  yields  a  more  substantial  reduction  of  the 
required  computational  effort  than  some  of  the  existing  acceleration  proce¬ 
dures. 

We  advocate  that  this  structure  should  be  analyzed  and  used  in  choosing  the 
appropriate  solution  procedure.  This  procedure  might  be  composed  by  blending 
several  of  the  acceleration  concepts  that  are  described  in  literature.  The 
four  test  problems  are  sketched  and  solved  with  several  successive  approxima¬ 
tion  methods.  These  methods  were  composed  after  analysing  the  structure  of  the 
problem.  The  required  computational  efforts  are  compared.  Among  other  things 
it  appears  that  combination  of  different  acceleration  techniques  is  sometimes 
counterproductive.  It  also  appears  that  it  often  is  important  to  select  the 
computational  method  on  the  basis  of  its  possibilities  for  efficient  handling, 
rather  than  on  its  convergence  rate. 
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On  Stationary  Reliability  Characteristics 
of  Complex  Systems  with  Repair 
by 

Peter  Franken 
Humboldt  University 
Berlin 

ABSTRACT 


1.  In  this  paper  a  unified  approach  is  given  to  derive  formulas  for  the 
state  probabilities, _  the  availability  and  the  interval  reliability  of  a 
complex  repairable  system  in  steady  state.  The  approach  is  based  on  an 
application  of  stochastic  processes  with  an  embedded  point  process,  cf. 
e.g.,  1 1  ]- [ 5 ] . 

Definition.  A  pair  [(X(t),-«=  <  t  <  °°) ,  <  n  <  °»)  ]  consisting  of  a 

stochastic  process  (X(t))with  the  state  space  (E,o(E))  and  a  sequence  of 
real  valued  r.v.'s  *(T  )  both  defined  on  the  same  probability  space 
(n,  Q)  and  satisfying  the  condition 


r  - ■— ->  ± 

n  n  -*■  ±” 


...  <  T„  <  0  <  T.  <  ...  Q-a.s. 


is  called  a  process  with  an  embedded  point  process  (abr.  PEP).  The  points 
T^  are  called  the  embedded  points. 

A  PEP[(X(t),(T  )]  is  called  synchronous  if  the  cycles  (X  (u)  =  X(T  +u), 
n  n  n 

0  _<  u  <  T^+^  -  T^),  <  n  <  ”,  form  a  stationary  sequence  (or  variables 

taking  values  in  a  suitable  space). 

A  PEP [ (X(t) ) , (T^) ]  is  called  stationary  if  for  every  u  c  R  the  shifted  PEP 

S  l(X(t)),(T  )]  =  [(X(t-u)),(T  -u)]  has  the  same  distribution  as  [(X(t),(T  )]. 
u  n  n  n 

We  use  the  notations 

dP  *  { P :  P  is  the  distribution  of  a  synchronous  PEP[ (X(t) , (I  ) )  with  the  property 
0  <  M P)  -  EpT1  .<  -} 

{ P :  P  is  the  distribution  of  a  stationary  PEP [ (X(t ) , (T  ) J  with  the  property 
0  <  A (P)  -  Ep[#(0,l)  n  (Tn))l  <  -} 


Theorem.  There  exists  a  one  to  one  mapping  between  the  families  &  and  &. 
For  a  given  P  c  &  Che  corresponding  P  is  determined  by 

P([(X(t)),(Tn)]  e  A)  =  A(P)'1  /  “  P({T1  >  u)  O  {Su<  [  (X(t) )  ,  (Tn)  ])  c  A)du  (1) 

From  (1)  we  obtain  the  so-called  stochastic  mean  value  theorem 
(f  _>  0  and  measurable) 

*  -1  T 

Epf (X(0) )  =  A(P)  *Ep</  lf(X(t))dt)  (2) 

This  formula  is  well-known  for  regenerative  processes,  i.e.,  if  the  cycles  of 

a  synchronous  process  are  i.i.d.,  cf  e.g.,  [6].  Formulas  (1)  and  (2)  are 

very  useful  in  queueing  theory,  cf  15). 

If  the  synchronous  PEP[ (X(t) , (T^) ]  is  metrically  transitiv,  one  can 

interpret  l(X(t)),(T  )]  and  the  corresponding  stationary  PEP | (X(t ) ) , (T  )] 

n  n 

as  two  equivalent  description  of  the  same  system  in  steady  state.  The  first 
one  corresponds  to  the  point  of  view  that  the  origin  is  an  arbitrarily  chosen 
embedded  point.  The  second  corresponds  to  the  origin  being  an  arbitrarily 
chosen  point  on  the  real  axis. 

2.  Let  us  consider  a  complex  system  with  repair.  We  assume  that  there  exist 
sets  G  and  B  (G-good  states,  B-bad  states)  and  that  CUB  =  E,  G  i  \  B  =  f'. 

The  system  is  failed  iff  its  state  belong  to  B. 

As  embedded  points  we  consider  the  instants  at  which  the  system  changes 
its  state  from  B  to  G.  If  the  system  is  in  steady  state  we  describe  its 
temporal  behavior  by  .a  synchronous  PEPl(X(t),(T  )]  and  the  corresponding 
(in  the  sense  of  Theorem  1)  stationary  PEP( (X(t ) ) , (T  ) ] . 

As  the  stationary  state  probabilities,  stationary  availability  and 
stationary  interval  relability  we  define  the  quantities 

Pj  -  P(X(0)«j),  A  -  P(X(0)  c  G),  Av  -  P(X(u)  e  G,  0  <  u  <  v),  (3) 

respectively.  The  lengths  of  generic  cycle  of  the  synchronous  process  can 


C18 


be  written  in  the  form  A  =  U  +  D,  where  Li  ,s  the  up-time  and  D  is  the  down¬ 
time  of  the  system  within  the  cycle.  Using  (1)  we  get 


v  EpU  +  EpD  v 


/  P(U  >  t )dt 


(4) 


EpU  +  EpD 


(5) 


Both  formulae  are  well  known  for  the  so  called  repairable  unit  if  the 
temporal  behavior  of  the  system  is  described  by  an  alternating  renewal 
process.  Formula  (3)  was.  proved  in  [7]  for  Markov  systems  with  finite 
state  space.  (We  point  out  that  the  formulas  (4),  (5),  as  well  as  the 
formulas  (6)  and  (7)  below  hold  for  an  arbitrary  state  space  E). 

3.  Let  C  e  o(E)  and  assume  that  the  instants  in  which  the  system  enters 
or  exits  the  set  C  belong  to  the  sequence  of  embedded  points.  Denote 


A  =  Ep(Tn+1  -  Tn),  *<L)  -  P (X(0  +  0)  c  L),  L  c  o(E>,  Ax  =  EpdJxfO  +  O)  =  x) , 

(we  assume  that  this  conditional  expectation  exists).  Using  the  formula 
(2)  for  these  embedded  points  and  f  =  lc  we  get 


1  /C  V(dx) 

P(X(0)  c  c)  =  -  /c  V(dx)  - 

E  x 


In  a  reliability  framework  (discrete  state  space)  we  have 

,  k  e  E 


"A 


L  it  A 
jeS  1  3 


(6) 


(6a) 


Consider  a  semi-Markov  process  with  an  arbitrary  state  space  E  such 
that  there  exists  a  unique  invariant  probability  measure  n  of  the  embedded 
Markov  chain.  If  we  consider  all  jump  points  as  embedded  point  the  formula 
(6a)  is  the  well-known  formula  for  the  stationary  distribution  of  a  semi- 


3 


Cl  9 

Markov  process,  cf  e.g.,  [1],  [2].  So  we  can  observe  that  the  formula  (5) 

holds  in  fact  for  arbitrary  PEP.  However,  unlike  semi-Markov  processes, 

the  determination  of  it  and  A^,  x  e  E  can  be  quite  difficult. 

We  get  another  formula  for  the  probability  P(X(0)  e  C)  if  we  consider 

the  entrance  points  into  C  as  embedded  points  (there  could  be  further 

embedded  points,  also).  We  denote  by  ((X(t),(Tn))  and  P  the  synchronous 

PEP  and  its  probability  distribution  and  use  the  notations  A^E-fT  -  T  ), 
r  •  1  p  n+1  n 

•  „  T 

ii  (L)  =P(X(0  +  0)  e  L),  L  e  o(E),  Ax  =  E^(/Q  1  lc  (X(t )  )dt  |  X  (0  +  0  =  x)  -  the 
conditional  mean  sojourn  -time  in  C  under  the  condition  that  the  process 
starts  in  x  e  C.  Using  (2)  we  get 

P(X(0)  e  C)  =  4  /_A*  it (dx)  (7) 

A  x 

For  discrete  E  we  consider  C  *  (k)  and  get 

p(k)  =  4  i(k)A^  (7a) 

The  formulas  (7)  and  (7a)  are  well-known  for  semi-regenerative  processes, 
cf  e.g.,  [1],  [2],  [4].  In  reliability  theory  one  can  use  formulas  (ba) 
and  (7a)  for  the  determination  of  the  availability  via 

A  -  E  p(k) , 
kcG 

cf  [1].  But  our  experience  shows  that  the  use  of  (5)  leads  to  much  simpler 
calculations. 
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Conservation  Laws  for  Queues 
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Daniel  P.  Heyman 
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Holmdel,  NJ 

ABSTRACT 


Queueing  systems  obey  a  conservation  law  which  relates  certain 
time  and  customer  averages.  This  talk  will  describe  a  generalization 
of  Little's  formula  (due  to  S.  Brumelle),  give  examples  of  its  use, 
and  sketch  a  new  proof. 


Queues  with  Non -Stationary  Input  Stream:  Rosses  Conjecture 

by 

Tomasz  Rolski 
Wrocjaw  University 

ABSTRACT 

Characteristics  of  queues  with  non-stationary  input  streams  are  difficult 
to  evaluate,  therefore  bounds  for  them  are  of  importance.  First  we  define 
what  we  understand  by  the  stationary  delay  and  find  out  stability  conditions 
of  single  server  queues  with  non-stationary  inputs.  For  this  purpose  we 
introduce  a  notion  of  ergodically  stable  sequence  of  random  variables.  The 
worked  out  theory  is  applied  to  single  server  queues  with  stationary  double 
stochastic  Poisson  arrivals.  Then  the  inter-arrival  times  do  not  form  a 
stationary  sequence  ("time  stationary"  do  not  imply  "customer  stationary"). 

We  show  that  the  averaqe  customer  delay  in  the  queue  is  greater  than  in  a 
standard  M/Gl  queue  with  the  same  average  input  rate  and  service  times.  This 
result  is  used  in  examples  which  show  that  the  assumption  of  stationarity 
of  the  input  point  process  is  non-essential. 
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